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Abstract 


Surface-dressed  optical  Bloch  equations  are  solved  for  a  two-level  atom  near 
or  adsorbed  on  a  rough  surface,  which  is  modeled  as  a  hemispheroid  protrusion  on 
a  perfectly-conduct ing  surface.  Effects  of  the  laser  bandwidth  are  included  by 
means  of  a  phase-diffusion  model.  The  presence  of  the  surface  roughness 
generates  a  dephasing  broadening  mechanism  on  the  adatomic  resonance 
fluorescence  spectrum,  which  can  be  strongly  enhanced  by  plasmon  resonances 
depending  on  the  shape  of  the  hemispheroid.  The  dependence  of  the 
electrodynamics  of  the  adatom  on  the  adatom-hemispheroid  distance  is  also 


evaluated. 


1 .  Introduction 


There  is  considerable  interest  in  the  interaction  of  laser  radiation  with 

adspecies  on  surfaces.*  In  particular,  researchers  have  studied  the  problem  of 

2-5 

surface-enhanced  spontaneous  emission  of  two-level  atoms  near  a  mirror  and 

decay  of  excited  molecules  near  a  small  metal  sphere.^  The  enhancement  of  Raman 
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scattering  and  fluorescence  by  small  metal  spheres  has  also  been  examined.  ’ 

However,  when  the  resonant,  driving  coherent  radiation  is  very  strong,  interesting 

9  10 

"resonance  fluorescence"  and  other  nonlinear  phenomena  can  occur.  *  In  these 

cases,  the  resonant  laser  puts  the  atom  or  molecule  in  an  environment  where  the 

probability  of  stimulated  emission  can  exceed  that  of  spontaneous  emission.  Under 

these  circumstances  we  find  back-transitions  and  Rabi  oscillations  of  state 

probability  amplitudes,  and  the  dynamic  Stark  splitting  of  resonances  and 

mutational  oscillations  of  the  emitted  light  intensity  become  important  parts  of 
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the  laser-driving  process.  *  In  a  strong  resonant  field,  nonlinear  scattering 

occurs,  and  multiphoton  processes  become  as  important  as  single-photon  processes, 

so  that  we  can  no  longer  rely  on  the  conventional  low-order  perturbation  theory  of 

fluorescence.  In  the  surface-free  case,  i.e.,  in  the  absence  of  a  solid  surface 

or  when  the  atom  is  very  far  from  the  surface,  resonance  excitation  of  a  two-level 

atom  by  a  laser  field  has  been  extensively  investigated,  using  the  powerful 
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optical  Bloch  equations.  '  Previously,  we  have  derived  a  set  of  surface-dressed 
optical  Bloch  equations* *  *"*  by  which  we  can  discuss  the  effects  of  the  surface- 
reflected  photons,  the  resonance  interaction  between  the  adatom  and  surface- 
plasmons,  the  collision  dephasing  of  the  adatom  produced  by  gas  atoms  in  the 
medium,  and  the  random  phase  fluctuation  of  the  (intense)  laser  field.  In  the 
present  paper,  we  extend  the  calculation  of  the  resonance  fluorescence  for  a  flat- 
surface^  ® 


to  the  case  of  an  interface  with  metallic  protrusions. 


Surface  roughness  can  be  characterized  as  a  set  of  protrusions  which  randomly 


(or  regularly,  in  the  case  of  a  grating)  stick  out  of  the  underlying  planar 
substrate.  While  the  protrusions  may  have  different  sizes  and  shapes,  in  this 
paper  we  model  the  protrusions  as  a  prolate  hemispheroid  on  top  of  a  plane, 
as  depicted  in  Fig.  1.  We  note  that  the  model  of  a  half  spheroid  on  a  flat 
perfect  conductor^  has  been  shown  to  be  identical  to  a  full  spheroid  in  a 
vacuum. Thus  the  following  discussion  can  be  also  used  for  ellipsoidal  cases. 
The  two-level  adatom,  which  is  located  at  a  distance  d  from  the  top  of  the 
hemispheroid,  is  driven  by  a  laser.  The  adatom  has  no  dipole  moment  in  its  ground 
state,  but  it  can  have  a  transition  dipole  connecting  the  ground  and  excited 
states.  The  emitted  photons  will  be  reflected  by  the  flat  mirror  and  the 
hemispheroid,  which  in  turn  provide  a  resonance  feedback  to  the  dynamics  of  the 
adatom.  The  spectrum  of  laser  light  scattering  in  the  presence  of  a  bump-plane 
surface  system  near  a  two-level  atom  is  then  the  subject  of  this  paper. 
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2.  Surface-Dressed  Optical  Bloch  Equations  for  a  Rough  Surface 


Let  us  evaluate  the  light  scattering  produced  when  a  laser  field  is  incident 

on  a  molecule  adsorbed  on  a  rough  surface.  The  surface  will  be  taken  to  be  a 

prolate  hemispheroid  protruding  from  a  flat  plane.  The  spheroid  is  assumed  to 

have  the  complex  dielectric  constant  e(w),  while  the  plane  is  taken  to  be  a 

perfect  conductor.  The  incident  laser  field  is  taken  to  be  propagating  along 

the  interface,  where  its  electric  field  is  taken  to  be  along  the  normal  of  the 
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interface.  The  prolated  spheroidal  coordinate  (£,  n,  $)  system  is  used  to 
calculate  the  reflected  field.  We  define 


and 


5,  -  (a  +  d)/f 
C0  -  a/f 


.  ,  2  .2.1/2 
f  -  (a  -  b  )  , 


where  a  and  b  are  the  semi-major  and  semi-minor  axes  of  the  hemispheroid 

(Fig.  1).  The  reflected  field  at  the  position  of  the  adatom  (dipole),  in  the 
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limit  of  the  near-field  approximation,  can  be  written  as 

Er  -  ~  ?  l'cn  W  +  - "  3 

r  n  n  n  1  Mf^)3 


(1) 


where  Q  denotes  the  Legendre  function  of  the  second  kind  and  u  is  the  induced 
dipole  moment  which  will  be  given  later.  The  expansion  coefficient  Cn  is  given 
by 


C 


n 


-  i)eo  :  *o  V. 
‘  eiQiuo> 


<ln  ♦  2 


(i  - 

"  EWW  -  Wn«0> 


(2) 


where  denotes  the  Legendre  function  of  the  first  kind  and  Eg  is  the  amplitude 
of  the  incident  laser.  The  induced  dipole  moment,  u,  of  the  adatom  is  given  by 

M  ’  aEtotal  * 

where  a  is  the  polarizability  along  the  semi-major  axis  and  Etotal  is  the  total 


electric  field  applied  to  the  adatom.  The  induced  dipole  moment  can  be 


expressed  by 


»-  “(-K  c„w 


3  +  E0^  * 


After  rearranging  Eq.(3)  and  substituting  back  to  Eq. (1),  we  obtain 

„  E0  ,  (1  -  ,  . , 

L  ~r\  t  r  \  __  r  r\%  t  r  \  ^  J  t 


r  1  -  r  £Q1(^q)  S0Qi ( ) 


(3) 


where 


1»(f51)' 


-  1)  l' 

t1  n 


<2n . 

EW  W  -  W  W 


As  shown  in  Eq.(4),  there  are  two  resonance  conditions:  1-r  approaches  zero 
and  e  equals  5qQ] (^1 )/Q1 (€Q) •  It  can  be  shown  easily  that  the  near-field 
approximation  will  break  down  when  the  factor  1-r  approaches  zero.  The  physical 
interpretation  of  this  factor  is  the  image  enhancement  effect.  The  importance 
of  the  image  effect  has  been  reviewed  by  Schatz.  Let  us  focus  on  the  second 
resonance  condition.  We  recall  that 

Q1  (O  =|  ln[|4-|]  -  1  . 

Hence,  we  can  write  the  second  resonance  condition  as 

U0/2)  m[(c0  +  1)/U0  -  1)]  -  5q/(5q-  1) 

E  =  U0/2)  ln[U0  +  1)/(50  -  1)]  -  1  •  (5) 


Using  the  identity 


u2  -  d'1  -l  C2n 


we  obtain 


e  =  1  - 


V?0)(V  ’> 


(6) 


In  the  limit  »  such  that  the  spheroid  tends  towards  a  sphere,  Eq.(6) 


becomes 


c(5q->“) 


3  C 


=  1  - 


^0  '  1 


-2 


(7) 


which  reduces  to  the  usual  surface  plasmon  resonance  condition  of  a  sphere, 
e  +  2  =  0.  We  will  discuss  below  how  to  use  this  reflected  field  to  calculate 
the  transition  dipole. 

To  treat  problems  of  excitation  and  dissipation  of  a  two-level  adatom  near  a 
hemispheroid,  a  self-consistent  approach  is  developed  in  which  the  dynamical 


behavior  associated  with  the  induced  transition  dipole  is  determined  by  surface 

11,12 


dressed  optical  Bloch  equations  (SBE).  ’  Essentially,  the  SBE  are  quantum 

A  A  A  A 

operator  equations  for  o^,  o21  and  W,  where  »  |  i>  <j  |  (i,j  -  1,2)  are  the 

A  A 

adatomic  transition  operators  and  W(t)  =  o22(t)  -  o^Ct)  is  the  population 
inversion  of  the  adatom.  By  defining 

(8) 
(9) 

where  is  the  laser  frequency,  we  can  write  the  rotating-wave  approximation 
(RWA)  for  the  SBE  of  the  adatom-bump  system  as 


S12(t)  «*  0-|2(t)  exp(ioiLt) 


S21(t)  =  o21(t)  exp(-iwLt) 


C  -s 

r 

V 

**  > 

s2,(t) 

-Y2  +  iA 

ifl”( t) 

2 

0 

i21(t) 

0 

a 

+ 

a 

dt 

W(t) 

a 

ia  (t) 

"Y1 

-iif  (t) 

W(t) 

"  Yi 

1 

0 

-R+( t) 

2 

-Y2  -  i A 

i12<t) 

0 

N. 

R±(t)  = 

R  exp[ + 

i 4>( t ) d  is  the 

time-dependent  Rabi  frequency, 

where 

fl  -  (2/ii)  l^il  ^(t)  and  E(t)  is  written  as  an  expectation  value  in  a  coherent 
state  of  the  laser  field  in  terms  of  the  phase  factor  <j>(t)  as 
E(t)  -  Eo^t)  expC-ioo^t  +  i <J>( t ) D  ♦  c.c.;  A  -  u>21  -  is  the  detuning;  the 
surface-induced  phase-decay  Yp  -  Yp  +  Y  ,  where  Y  is  determined  by  the 


reflected  field  which  is  emitted  by  the  induced  adatom  and  reflected  by  the 

11  12 

bump-plane  interface.  In  the  case  of  a  rough  surface,  Yg  is  given  by  * 

Y  =(2/tf)lm(F)  ,  (11) 

o 

where 

(1  -  e)  CnQ{  U, ) 

F '  '_:T  C£  VV  -  e,a)(«0)  *  r  ]  • 

To  include  the  effects  of  the  laser  bandwidth  Y^,  the  phase-diffusion  model 
for  the  laser  field  has  been  used  by  defining  the  following  correlation  function 
for  the  Rabi  frequency: 

<<  Q  (tj)  S1+(t2)  >>  =  Q2  exp(-YL  |t2  -  tj  | )  ,  (12) 

where  the  double  bracket  signifies  two  averages:  one  is  over  the  stochastic 
ensemble  and  the  other  is  a  quantum  mechanical  average. 

a 

In  the  weak-field  or  large-detuning  cases,  where  W(t)  -  <<  W(t)  >>  *  -1,  the 
population  inversion  of  the  two-level  adatom  has  the  following  analytical  form 
(for  t£  0) : 

Y+Y  Y-Y-Y 

w(t)  -  -i  *  a2  1 - ^ - 1 - j - 2l  2  exp( -Y.  t) 

Y1[(Y2  +  V  *  A  3  V(Y1  '  Y2  "  V  +  A  3 

[(Y2  +  Yl)(Y1  -  Y2  -  Yl)  +  A2]  cos(At)  +  A(2Y2  -  2Yl  -  Y] )  Sin(At) 
C(Y2  +  YL)2  +  a2][(7i  ~  \  ~  \)Z  + 

x  exp[-(Y2  +  YL)t]  }  .  (13) 

The  corresponding  power  spectrum  S(u)  for  the  scattered  light  of  the  laser  can 
be  calculated  through  the  Fourier  transform  of  the  dipole-dipole  correlation 

A  A 

function  <<S21 (t2)S12(t1 )>>: 


where  the  parameter  a  is  defined  as 


Y2/2  {2 


(16) 


In  Eq . (15),  we  have  assumed  zero  laser  bandwidth,  Y.  =0.  The  detuning 


parameter  is  defined  as 
52'  *=  A  ♦  A 

c 

where 


(17) 


as  =  Ain  ♦  (j)/a)2]1/2  -  i } 

and  the  widths  of  the  profiles  are 

2  2 

S  =  2Y  (— _ *  —  ) 

0  2  1  rt2  2 ' 

52  +  A 


(18) 


(19) 


S±  =  s+  =  s_  =  y2C 


(1  +  q)52  2aA 

2  2 
52  +  A 


(20) 


Eq.(15)  was  derived  under  the  assumption  of  a  very  strong  coherent  laser  field, 
2  2-2 

i.e.,  52  +  A  >>  Y2.  The  strong-field  spectrum,  Eq.(15),  has  three  spectral 

peaks:  the  coherent  and  incoherent  parts  of  the  Rayleigh  scattering  centered  at 
a)  =  oiL  and  the  two  side  bands  centered  at  u  *  +  52'  due  to  a  multiphoton 

process  and  fluorescece  of  the  adatom,  referred  to  as  the  three-photon 
fluorescence  components,  respectively.^  We  shall  see  in  the  next  section  how 
the  resonance  coupling  between  adatom  and  bump  influences  the  strong-field 
three-peak  spectrum. 


3.  Numerical  Results 

Using  Eqs.(13),  (I1!)  and  (15),  which  describe  the  excitation  and  relaxation 
behavior  of  a  two-level  atom  near  or  adsorbed  on  a  metal  surface  with  roughness 
(a  bump  on  a  flat  surface),  we  have  evaluated  numerically  some  typical  cases  for 
the  rough  surface-modified  resonance  fluorescence  spectrum.  In  all  the 
calculations,  the  laser  photon  energy  is  taken  to  be  2.75  eV. 


Figure  2  displays  the  rough  surface-induced  phase-decay  constant  Y  as  a 

s 

function  of  the  semi-minor  axis  b  of  the  surface  protrusion,  with  the  semi-major 
axis  a  for  the  surface  protrusion  fixed  at  100  A.  There  is  a  sharp  resonance 
peak  at  b  =  38  A,  corresponding  to  a  resonant  excitation  of  plasmons  in  the  bump 
which  can  strongly  enhance  adatom-bump  interaction,  through  the  reflected  field 
at  the  atomic  site.  Figures  3(a)  and  (b)  display  the  influence  of  the  plasmon 
resonance  in  the  bump  on  the  adatomic  population  distribution.  The  time 
oscillations  of  the  population  will  decrease  as  the  semi-minor  axis  b  of  the 
spheroidal  bump  approaches  the  plasmon  resonance  b  =  38  A  [Fig.  3(b)], 
indicating  that  the  enhancement  of  the  adatom-bump  coupling  leads  to  more 
damping  of  the  adatomic  oscillator. 

Figure  4  displays  the  resonance  fluorescence  spectrum  in  the  weak-field 
driving  case  [Eq.(lH)],  consisting  of  the  Rayleigh  and  fluorescence  peaks.  The 
former  peak  in  our  model  is  determined  by  the  laser  bandwidth  Y^.  The  latter 
peak  (fluorescence)  strongly  depends  on  the  adatom-bump  distance  and  the 
resonance  coupling  between  the  adatom  and  the  bump.  When  the  detuning  A  gets 
smaller  and  the  semi-minor  axis  b  of  the  bump  spheroid  reaches  the  plasmon 
resonance  condition  [see  Fig.  Mb)],  the  two  peaks  will  overlap,  where  the 
fluorescence  peak  is  strongly  broadened  by  resonant  adatom-bump  coupling.  In 
this  case,  the  height  of  the  Rayleigh  peak  will  depend  on  the  atom-bump  distance 
and  parameter  b. 

Figure  5  displays  the  typical  three-peak  spectrum  of  the  strong-field-driven 
adatom  using  the  formula  of  Eq .(15) .  In  Fig.  5(a)  we  can  distinguish  the 
Rayleigh  (central)  peak,  the  three-photon  peak  (left)  and  the  fluorescence  peak 
(right).  In  these  pictures  we  have  drawn  just  the  profile  of  the  incoherent 
component  of  the  Rayleigh  scattering,  as  the  coherent  component  will  not  be 
broadened  by  the  surface  bump. 
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In  the  smaller  detuning  case,  i.e.,  Fig.  5(a),  the  spectrum  has  a 
distinctive  three-peak  nature.  In  this  case,  the  three-photon  peak  (left)  has 
measurable  height,  which  is  almost  comparable  with  the  fluorescence  peak.  It  is 
seen  that  the  smaller  the  adatom-bump  distance,  the  larger  the  broadening 
equally  in  the  three  spectral  peaks,  due  to  the  enhancement  of  the  adatom-bump 
coupling  for  smaller  distances.  In  the  larger  detuning  case,  i.e.,  Fig.  5(b), 
the  driving  effect  of  the  applied  laser  field  is  weaker  due  to  the  large 
detuning  A.  The  height  of  the  three-photon  peak  is  decreased  and  the  three-peak 
spectrum  is  transformed  to  the  weak-field  two-peak  structure.  In  this  case, 
smaller  adatom-surface  distances  still  lead  to  larger  broadening,  as  in  the 
former  case. 

Summary 

A  simple  model  has  been  developed  to  describe  the  influence  of  surface 
roughness  on  the  dynamic  behavior  of  a  two-level  adatom  near  a  metal  surface. 
Surface-dressed  optical  Bloch  equations  (SBE),  which  account  for  surface- 
reflected  photons  and  the  resonance  interaction  due  to  plasmon  excitation  in  the 
metal  protrusion,  have  been  derived.  Based  on  the  SBE,  the  population 
distribution  of  the  adatom  and  the  power  spectrum  of  the  scattered  light,  as 
well  as  the  surface-enhanced  dephasing  broadening  of  the  atomic  spectrum,  have 
been  evaluated.  Collisions  with  foreign  gas  atoms  and  effects  of  the  laser 
bandwidth  can  also  be  considered  in  our  model.  It  is  shown  that  the  spectrum  of 
scattered  light,  as  well  as  the  level  population  of  the  adatom,  are  strongly 
enhanced  by  the  resonantly  excitation  of  plasmons  in  the  hemispheroid  on  the 
metal  surface.  The  displayed  surface  enhancement  for  the  adatom-bump  system 
should  be  very  useful  in  understanding  the  spectroscopy  and  chemistry  associated 


with  rough  surface. 
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Geometry  of  the  surface  protusion.  The  semi-major  axis  is  a  and  the  semi¬ 
minor  axis  is  b.  The  prolate  spheroidal  coordinates  are  defined  by  a  scale 

2  2  1/2 

factor,  f,  given  by  f  *  (a  -  b  ) 

Rough  surface  induced  phase-decay  constant  Yg  as  a  function  of  the  semi¬ 
minor  axis  b  of  the  surface  protusion,  where  the  semi-major  axis  is  fixed 
at  100  A.  Ys  is  in  the  unit  of  Einstein's  spontaneous  decay  constant  A. 
Time  evolution  of  the  population  inversion  in  the  weak-field  or  large- 
detuning  limit  with  (Y^,  0,  A)  =  (0.3,  0.05,  5)  in  the  unit  of  Einstein's 
spontaneous  decay  constant  A  for  b  =  34  A  [Fig.  3(a)]  and  b  -  38  A 

[Fig.  3(b)],  w;tK  a  fixed  at  100  A.  The  horizontal  axis  corresponds  to  time 
-1  2 

in  the  unit  A  .  The  quantity  [w(t)  +  1  ] / SI  is  shown  along  the  vertical 
axis.  The  curves  numbered  1 ,  2,  3  and  4  correspond  to  the  reduced  atom- 
bump  distance  d  =  20,  25,  30  and  35,  respectively  ( d  »  2nd/A,  where  A  is 
the  laser  wavelength;. 

Resonance  fluorescence  power  spectrum  of  scattered  light  with  (Y^,  Q,  A)  - 
(0.3,  0.05,  5)  in  the  unit  A  for  b  =  40  A  [Fig.  4(a)]  and  b  «  38  A 
[Fig.  4(b)].  Curves  1,  2  and  3  correspond  to  the  reduced  atom-bump 
distances  d  =  20,  25  and  30,  respectively.  X  denotes  (w  -  u)L)/(u)21  -  w^) 
on  the  horizontal  axis.  The  units  along  the  vertical  axis  are  arbitrary. 
Strong-field  surface-modified  resonance  fluorescence  spectrum,  where  just 
the  incoherent  component  of  Rayleigh  scattering  is  included,  with 

(A,  0,  b)  =  (1,  10,  38)  for  Fig.  5(a)  and  (A,  0,  b)  -  (10,  10,  38)  for  Fig. 

Arc  ir\  And 

5(b),  A  and  0  ^  the  unit  of  A f> b  in  the  unit  A.  Curves  1,  2  and  3 

correspond  to  the  reduced  atom-bump  distances  d  -  20,  25  and  30, 
respectively.  X  denotes  (w  -  on  the  horizontal  axis.  The  units 

along  the  vertical  axis  are  arbitrary. 
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